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We investigate theoretically and numerically the photonic band structure in the optical domain 
of an array of core-shell metal-semiconductor nanowires. Corresponding negative-index photonic 
bands are calculated, showing isotropic equifrequency surfaces. The effective (negative) electric 
permittivity and magnetic permeability, retrieved from S-parameters, are used to compare the per¬ 
formance of such nanowire arrays with homogeneous media in canonical examples, such as refraction 
through a prism and flat-lens focusing. Very good agreement is found, indeed confirming the effec¬ 
tive medium behavior of the nanowire array as a low-loss, isotropic (2D) and bulk, optical negative 
index metamaterial. Indeed, disorder is introduced to further stress its robustness. 

PACS numbers: 78.20.Ci,78.67.Pt,42.25.Bs,42.70.Qs 


I. INTRODUCTION 

Experiments in composite systems looking for a neg¬ 
ative index of refraction (both permittivity and perme¬ 
ability being negative) at gigahertz frequencies were the 
first of this kind aiming at controlling the flow of elec¬ 
tromagnetic waves; indeed, some fascinating properties 
have been since explored based on such compositeS^. 
However, to scale down these properties, not found in 
natural media, from the gigahertz regime to the opti¬ 
cal domain, has posed a major challengeespe cially 
with regard to the negative magnetic permeability 5 16i 
(negative electric response is trivial in the optical do¬ 
main with metals). Despite some drawbacks of the orig¬ 
inal design, as for instance anisotropy, the initial philos¬ 
ophy to achieve this goal was to miniaturize the original 
ones, namely, split-ring resonators, magneto-dielectric el¬ 
ements, and related 10 12 13 . Nonetheless, the saturation 
of the magnetic response at high frequencies restricts this 
approach to the far infrared^. Other desig ns hav e been 
proposed that somehow mimic their behavio j 12 * 17 ! ®, one 
of the most widely employed being metallic fishnets 21 . 
The emergence of ohmic losses in metals in the optical 
domain hinders also their performance. 

To overcome this problem, various attempts have 
been made by exploiting magnetic resonances occurring 
in structures made of high-permittivity materials®®, 
some of them combined with secondary plasmonic struc¬ 
tures providing the electrical response 32 3 A Very re¬ 
cently, semiconductor nanostructures exhibiting mag¬ 
netic resonances have been proposed as possible non¬ 
absorbing metamaterial elements as well as emitters 
with unusual properties 36 ®. Actually, simple arrange¬ 
ments based on this idea have been proposed such as 
core-shell 3D and 2D metallo-dielectric nanostructures, 
respectively, nanospheres and nanowires, geometrically 
tuned so that their two first order Mie resonances (mag¬ 


netic and electric) spectrally overlap, which might play 
the role of doubly-resonant meta-atoms®!®_directional 
scattere d 44 * 45 *, and transparent nanowires 4 ^®. Be¬ 
side their robustness and isotropy stemming from their 
highly symmetric, single-structure meta-atoms, a the¬ 
oretical figure of merit ( f.o.m . = — 5R(n e ff)/^(n e ff)) of 
the order of f.o.m. ~ 200, has been rep orted in the 
case of core-shell nanowire (NW) arrays®. Incidentally, 
NW arrays have been proposed in turn as hyperbolic 
metamateriald®. 

Nevertheless, a full description of their wave propa¬ 
gation properties must be done to determine whether 
their negative refraction behavior is a photonic-crystal 
diffraction effect that depends on the geometry sam¬ 
ple, or a fully isotropic effective negative-index behavior 
(which we will referred to as negative-index metamate¬ 
rial, NIM). In fact, it has been pointed out in Ref. f50l that 
the constitutive parameters retrieved through effective 
medium theories in other dielectric-rod ordered arrays 
do not reproduce the negative refraction observations in 
those composites, being mostly due to photonic crystal 
diffraction®*®. In this work, we thoroughly investigate 
the photonic band structure of an array of core-shell 
metal-semiconductor NWs, which was shown in Ref® 
through parameter retrieval to exhibit low-loss, negative 
index behavior. This allows us to rigorously demonstrate 
(through equifrequency surfaces) and improve (fine tun¬ 
ing core-shell geometrical parameters) its highly isotropic 
character. Indeed, taking advantage of high isotropy and 
low losses, flat-lens focusing is presented through thick 
slabs with subwavelength resolution and dipole position 
independence, in good agreement with the expected be¬ 
havior for a homogeneous slab with effective parameters 
as retrieved for the NW array. Limitations are also inves¬ 
tigated with regard to image resolution, photonic-crystal 
surface termination, and disorder. 

This paper is organized as follows. Corresponding 
negative-index photonic bands and equifrequency sur- 
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FIG. 1. (Color online) Schematic of the photonic crystal 
configuration under study, an hexagonal array of core-shell 
nanowires, including the relevant geometrical parameters. 
Only propagation in the plane perpendicular to the (infinitely 
long) NW axis will be considered with transverse electric po¬ 
larization. 

faces are calculated in Sec. II, confirming the effective 
medium behavior. In Sec. Ill, the performance of such 
NW arrays in the canonical example of flat-lens focusing 
through a thick slab, as compared with a homogeneous 
slab with the retrieved electric permittivity and magnetic 
permeability. This is also done for plane wave refraction 
through a prism in Sec. IV, introducing in turn disorder 
to support the robustness of the NW array as a low-loss 
optical NIM. Finally, Section V contains our concluding 
remarks. 


II. PHOTONIC BAND STRUCTURE 

First of all, we characterize the electromagnetic wave 
propagation in the 2D periodic composite medium con¬ 
sisting of core-shell metal-dielectric (infinitely long) NWs 
in a regular hexagonal array (see Fig. [lj. For this pur¬ 
pose, we plot in Fig. [2] the 2D photonic band structure 
for the maximum symmetry directions for TE-polarized 
waves (magnetic field along the NW axis) of an hexago¬ 
nal array with a lattice parameter of a = 350 nm made 
of Ag@Si core-shell NWs with external radius R 0 ut = 
170 nm, and inner radii R[ n = 80,90 nm. Details on the 
FEM numerical calculations of the photonic band struc¬ 
ture are given in the appendix. Also included in Fig. [2] 
is the total Mie scattering efficiency of a single Ag(Q)Si 
NW, obtained from Mie formulae for infinitely long core¬ 
shell cylinderd^, explicitly plotting the contributions of 
the first three terms (Qo, Qi, Q2) stemming, respectively, 
from the zero-order magnetic, first-order (dipolar) elec¬ 
tric, and second-order (quadrupolar) electric terms in the 
Mie multipolar expansion. The first one is due to a strong 
circulation of the electric displacement field inside the di¬ 
electric shell, so the position of the resonance mostly de¬ 
pends on the shell external radius. The Q 1 contribution 
corresponds to the localized surface plasmon resonance of 
the metal core. For obvious reasons, i.e. NIM behavior 4 -^, 
we choose the appropriate geometrical parameters so that 
both (electric and magnetic) resonances spectrally over¬ 
lap, as clearly seen in Fig. [2] indeed, another inner radius 
Ri n = 90 nm has been considered that slightly improves 


FIG. 2. (Color online) In-plane photonic band structure (cen¬ 
ter) of an hexagonal array (a =350 nm) of Ag@Si core-shell 
nanowires with external radius Rout = 170 nm and internal 
radii: (top) Ri n = 80 nm and (bottom) Ri n = 90 nm. In 
both plots: (left) The Mie scattering efficiency for an isolated 
nanowire is also shown (black curve), including separately 
the contributions from the first three multipolar terms (red, 
green, and blue curves); (right) phase index n p corresponding 
to the first two bands along the main TM (red curves) and 
TK (blue curves) directions. 

such NIM behavior, as will be shown below. In addition, 
the lattice constant is such that the high filling fraction 
condition needed to achieve doubly negative effective pa¬ 
rameters is fulfilled^ the lattice choice, hexagonal, fa¬ 
vors this close packing and in turn isotropy. 

The photonic bands present a partial gap in the TM di¬ 
rection in the first band and a full gap between the third 
and fourth bands. Besides, the slope of the second one is 
negative, and corresponds to a mode propagation with a 
negative phase velocity, as evidenced in the right graphs 
in Fig.[2j where the phase index n p corresponding to the 
first two bands along the main TM and TK directions 
is plotted. Furthermore, this band lies on the spectral 
region where the magnetic and electric resonances over¬ 
lap, making it specially relevant to elucidate the effective 
NIM behavior stemming from strong individual electro¬ 
magnetic responses of the isolated NWs. Thus let us 
study in detail this negative-phase second band in Fig. [2j 

To this end, we represent in Fig. j3^a) the lowest 
four photonic bands, upon superimposing the curves at 
various propagation directions in between the maximum 
symmetry directions, TM and TK (red and blue curves, 
respectively), covering all relevant directions within the 
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FIG. 3. (Color online) (a) Photonic band structure as in 
Fig. [5] (a = 350 nm, R ou t = 170 nm), but superimposing vari¬ 
ous (30) bands at different directions within the first Brillouin 
zone in between TM (red curves) and TK (blue curves), (b) 
Equifrequency surfaces for several frequencies (dashed hori¬ 
zontal lines in (a)) covering mainly the second band. In both 
(a,b): (top) Rin = 80 nm and (bottom) R in = 90 nm. 
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FIG. 4. (Color online) Spectral dependence (within the sec¬ 
ond NIM band) of the refractive index n e ff (a), impedance Z e & 
(b), effective dielectric permittivity £ e ff (c) and magnetic per¬ 
meability /x e ff (d) retrieved from the S-parameters for Ag@Si 
core-shell nanowires, arranged in a hexagonal lattice with lat¬ 
tice parameter a — 350 nm, with outer radius R out = 170 nm, 
and two inner radii: Ri n — 80 nm (blue lines) and Ri n = 90 
nm (red lines). Solid curves: real parts; dashed curves: imag¬ 
inary parts. The phase index n p obtained from the photonic 
bands in Fig. [5] is also plotted (with symbols, top left). 


first Brillouin zone. In this manner, the partial gaps and 
the full gap at 300 THz are clearly observed. More impor¬ 
tantly, the nearly isotropic character of the second band 
is inferred, improving upon slightly increasing the NW 
inner radius from a dispersion of 10 THz (^ 5%) with 
wavevector direction shift for R in = 80 nm, to a mere 2 
THz (~ 1%) for R in = 90 nm. 

In order to clarify isotropy of the wave propagation 
properties in the second band we plot in Fig. §b) its 
equifrequency surfaces (EFSs) for several frequencies 
[horizontal lines in Fig. |3](a)] , bearing in mind the hexag¬ 
onal symmetry of the crystal. In the beginning of the sec¬ 
ond band, at 215 THz (1.4 /im), we have the well known 
pockets due to the partial gap along the TM direction, 
whereas the only allowed propagation modes are close to 
the TK directions. For frequencies above the partial gap 
[see Fig. [3|a)], all the mode propagation directions are 
again available and the EFSs become more circular, as 
shown at 235 THz (1.28/im): Isotropy is indeed remark¬ 
able for the NW R[ n = 90 nm inner radius at 245 THz 
(1.22 /im), with a nearly circular EFS. Then, as the fre¬ 
quency is further increased, at 270 THz (1.11 /im) and at 
280 THz (1.07 /im), the EFSs of the second bands shrink 
and deform (become more hexagonal), overlapping with 
the third-band pockets appearing along the TK direc¬ 
tions, thus breaking isotropy. Overall, it should be em¬ 
phasized that there is a wide frequency range (215-260 
THz) in which the EFSs are almost circular (specially for 
Ri n = 90 nm) and the crystal might be described as an 
isotropic negative-index medium. 

In the light of these results, it seems possible to make a 


fully isotropic 2D NW photonic crystal with negative in¬ 
dex of refraction by taking advantage of both the negative 
index of the second band, due to the strong electromag¬ 
netic resonances of the isolated NWs, and of the isotropy 
that the photonic NW arrays shows within such band. To 
this aim, we first show in Fig.[4]the retrieved effective pa¬ 
rameters (S-parameter retrieval procedur e] 55 * 56 *, through 
FEM calculations as detailed in the appendix) for the 
electric permittivity and magnetic permeability, along 
with the resulting index of refraction and impedance: 
The doubly negative behavior with $ft(£ e ff), 5ft(/i e ff) < 0 
within the NIM band is evident in Fig. [4j Low losses 
are also obtained: 9(£ e ff), 9(/i e ff) ~ 0 [note that both 
might be slightly negative, in principle considered un¬ 
physical, but in practice as a result of the limitations 
of the parameter retrieval procedure with 9( n e ff) ~ 0]. 
The effective index is negative within the NIM band with 
nearly negligible losses: Sft(n e ff ) < 0, ^(n e ff) ~ 0. Indeed, 
the retrieved refractive index for R[ n = 90 nm yields an 
extremely large value of the f.o.m. ~ 200. Recall that 
the absorptive losses of a single core-shell NW are very 
smalP^, so that the remaining (small, but non-negligible) 
NIM losses are very likely due to scattering. The phase 
index n p obtained from the photonic bands in Fig. [2] is 
also plotted, in good agreement with Sft(n e ff). 

Incidentally, note that s e R and /i e ff do not become -1 at 
the same frequency, which results in an impedance Z e ^ = 
(/igff/^eff) 1 / 2 > 1. Further tuning of the core-shell NW 
parameters could be done in order to achieve the ideal 
Z e ff = 1; despite that, the resulting impedance is only 
slightly larger than 1 so that reflections at the interface 
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FIG. 5. (Color online) Transmission of the electromagnetic 
fields produced by an in-plane point dipole source at 235 THz 
(1.28 /im) located at 1.8 /xm above a slab made of Ag@Si 
core-shell nanowires with inner and outer radius Ri n — 90 
nm and Rout = 170 nm, respectively, arranged in a hexago¬ 
nal lattice with lattice parameter a = 350 nm: Slab surface 
along the TM direction. Left: Norm of the in-plane electric 
field for a dipole parallel (top, p x ) and perpendicular (bot¬ 
tom, p y ) to the surface slab respectively. Right: Norm of the 
Poynting vector for a dipole parallel (top) and perpendicular 
(bottom) to the surface slab, respectively. The results for a 
homogeneous slab are also shown, with £ e ff = — 0.734 — zO. 157 
and jjLefE = —2.094 + z0.532, so that n e & — —1.264 + x0.0069 
(f.o.m. ~ 180) as obtained in Fig. [4] at 1.28 /xm. 


are small, as will be shown below. However, at this point 
the important issue needs to be addressed whether or 
not this structure behaves as an isotropic homogeneous 
medium. 


III. FLAT-LENS FOCUSING: ISOTROPIC NIM 
BEHAVIOR 


To verify that the core-shell NW array presents an 
isotropic negative index of refraction, rather than con¬ 
sidering plane wave negative refraction (which has been 
confirmed up to large angles of incidence in Ref. [42]) , we 
study the imaging of a point source through a slab, which 
no doubt places a more stringent test for it involves all 
plane wave angular components. We carry out numeri¬ 
cal simulations of the transmission of the electromagnetic 
field radiated by an in-plane point dipole source at 235 
THz (1.28 fi m) placed in front of a finite (22 NW thick, 
namely, 6.7 /xm long) slab made of the structure studied 
in the previous section (see the appendix for the param¬ 
eters of the FEM calculations); we focus on the R[ n = 90 
nm case, which exhibits slightly higher isotropy. In Fig. [5] 
we plot the norm of the electric field (left column) and 
the Poynting vector (right column), respectively, for a 
dipole parallel to the surface of the slab (top), clearly re¬ 
vealing the effects of a medium that possesses a negative 
index of refraction. At the slab interfaces, light undergoes 
negative refraction and its field components are focused 
twice, one inside and other outside of the slab: The image 
produced on the other side of the slab presents similari¬ 
ties with the real one. For a dipole perpendicular to the 
slab (see Fig. [51 bottom), again, we see the same phe¬ 
nomenology of negative refraction through the slab with 
the expected focusing effect for this dipole orientation. 
In addition, it is remarkable the fact that the Poynting 
vector is always zero along the propagation direction par¬ 
allel to the dipole moment axis, so that light undergoes 
no significant spurious scattering losses when traveling 
through the photonic crystal; which in turn contributes 
to the effective medium parameters as losses, reportedly 
extremely small in this NIlVpSl, and thus in agreement 
with the above simulations. 

There fore the main features of flat lensing are fully 
retrieved 3 5 * 6 * 14 l Furthermore, the agreement is remark¬ 
able with the results (shown also in Fig. [5]) obtained upon 
considering an isotropic and homogeneous slab with ef¬ 
fective parameters at 1.28 /xm given in Fig. [4j it should 
be emphasized that the (linear) color scale of both the 
electric field (left) and Poynting vector (right) norms are 
identical in all plots (columns), either for the NW array 
or for the homogeneous medium, and also that the scale 
is neither saturated nor altered within the slabs. Re¬ 
markably, note that very similar transmitted foci are pro¬ 
duced: this is explicitly shown in Fig. [6] where the power 
transmitted at the image plane is presented for both 
the NW array and the homogeneous lenses. The full- 
width at half-maximum is of the order of 0.65 /xm^ A/2, 
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FIG. 6. (Color online) Poynting vector norms from Fig. [5] at 
the image plane. 


yielding subwavelength resolution. Of course, since the 
source is not located deep in the near-field, superlensing 
in the proper sensed cannot be observed and resolution is 
diffraction limited (ultimately, resolution should be lim¬ 
ited for the NW array to NW diameter). Overall, we 
clearly observe the phenomenology of (all-angle) nega¬ 
tive refraction at the interface of the slab and the flat- 
lens focusing process. This supports the highly isotropic, 
low-loss NIM behavior of the core-shell NW array, and 
the accuracy of the retrieved parameters, in turn ruling 
out photonic crystal collimation effects^®. 

Nonetheless, we plot in Fig. [7] the EFSs for the partic¬ 
ular case shown in Fig.[5j Momentum conservation along 
the impinging surface (wavevector component along TM) 
is explicitly denoted by the green dashed line for a par¬ 
ticular angle of incidence. Such plane wave component 
should be properly refracted with negative phase velocity 
(blue vector) and positive Poynting (red) vector. Since 
the EFS is nearly circular and larger than the air light 
cone, even plane wave components impinging from vac¬ 
uum at grazing angles are negatively refracted. Note also 
that no additional diffracted components are expected. 
Therefore, all components of the dipole source should 
undergo negative refraction in a consistent manner; this 
corroborates the above presented dipole focusing through 
a slab (see Fig. [5|. 

Such NIM focusing behavior is in turn preserved upon 
moving the source. First, we show in Fig. [8] the negli¬ 
gible impact of shifting the dipole source (either p x or 
p y ) by half-period parallel to the surface, as compared 
to Fig. [5] and in agreement with the homogeneous slab, 
further ruling out any artifact induced by single meta¬ 
atoms. Likewise, dipole positions are shifted perpendic¬ 
ular to the surface in Fig. [9| either closer (top) or farther 
(bottom); since this source displacement is more involved 
than the parallel one above, results are also included for 



FIG. 7. (Color online) EFS at 235 THz (1.28 pm) (blue cir¬ 
cle) of a hexagonal lattice (lattice parameter a — 350 nm, 
Brillouin zones delimited by red dashed hexagons) of Ag@Si 
core-shell nanowire with inner and outer radius R lu = 90 nm 
and Ro Ut = 170 nm. EFS of light impinging from vacuum/air 
is shown as a black circle. A dashed green (vertical) line de¬ 
notes wavevector conservation along the TM direction for a 
particular angle of incidence, refracted accordingly as indi¬ 
cated by phase/energy (blue/red, respectively) vectors. 


the homogeneous slab. According to the canonical NIM 
focusing effect (keeping the slab thickness fixed), if the 
dipole source is closer (respectively, farther) to the slab, 
the virtual image inside the slab gets closer (respectively, 
farther), whereas the image behind the slab appears far¬ 
ther (respectively, closer) from the slab surface; this is 
properly reproduced by the numerical calculations for the 
Ag@Si core-shell nanowire array. 

Nonetheless, we would like to stress that the direction 
along which the array is terminated at the surface plays 
a crucial role despite isotropy. We calculate in Fig. [10] 
the dipole source transmission through a slab identical 
to that used in Fig. [5j except for the fact that the slab 
surface is along the TK direction: it is evident that NIM 
focusing is remarkably worse. Two effects contribute to 
hinder it. First, surface roughness is larger along the TK 
direction, thus leading to larger scattering losses from sin¬ 
gle surface nanowires, as observed in Fig. [lOf through both 
larger electromagnetic fields on the slab surface close to 
the source, and larger reflection losses. Second, momen¬ 
tum conservation parallel to the TK direction leads to a 
scenario different from that along TM: apart from the 
expected negatively refracted beam, another refracted 
beam appears beyond a relatively small incidence angle 
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FIG. 8. (Color online) Top: Transmission through a slab 
made of a hexagonal array of Ag@Si core-shell nanowires (a = 
350 nm, Ri n = 90 nm, Rout = 170 nm, TM surface direction), 
of the electromagnetic fields produced by an in-plane point 
dipole source at 235 THz (1.28 /mi), as in Fig. [5] (Poynting 
vector only), but for vertical/horizontal dipoles shifted half¬ 
period parallel to the slab surface as depicted (bottom). The 
results for the equivalent homogeneous slab are identical to 
those in Fig. [5] and thus not shown. 


as a diffracted beam at nearest Brillouin zone. This im¬ 
plies that dipole source components at a certain angle 
with respect to the surface normal exhibit two diffracted 
beams that preclude proper NIM focusing: such behavior 
is evidenced in Fig.[l0]by crossing beams inside the slab. 


IV. REFRACTION THROUGH A PRISM: 

TOLERANCE TO DISORDER 

In the preceding section we have demonstrated by 
means of flat-lens focusing that this doubly-resonant NW 
array behaves effectively as an actual isotropic (bulk) 
optical NIM when arranged in a close hexagonal lattice 
packet, stemming from the isotropic negative-phase be¬ 
havior of the second photonic band. However, it has been 
shown that ordered configurations that exhibit a nega¬ 
tive index of refraction may lose such effective property 
when the geometry is disordered. In this regard, it should 
be recalled that it has been shown in other cylinder ar¬ 
rays that Mie resonances (rather than Bragg scattering) 
govern in part the appearance of spectrally over lapping, 
photonic band gaps and higher-order bands^^. There¬ 
fore, if Mie (rather than Bragg) scattering underlies the 
negative-index behavior, the latter should be more toler- 


v = 235 THz A = 1.28 nm 
Dipole p x Dipole p y 



FIG. 9. (Color online) Transmission through a slab made of a 
hexagonal array of Ag@Si core-shell nanowires (a = 350 nm, 
Rin = 90 nm, R out = 170 nm, TM surface direction), of the 
electromagnetic fields produced by an in-plane point dipole 
source at 235 THz (1.28 /mi), as in Fig. [5] (Poynting vec¬ 
tor only), but for horizontal (left column) and vertical (right 
column) dipoles shifted perpendicularly to the slab surface, 
both down (top) and up (bottom), as depicted by the left 
schemes. Results are also included for a homogeneous slab 
with e e ff = —0.734 — z0.157 and /x e ff = —2.094 + z0.532, so 
that n e ff = —1.264 + z0.0069 ( f.o.m . ^ 180) as obtained in 
Fig- 0 at 1.28/xm. 


ant to disorder than usual photonic crystals. 

In order to show the robustness of doubly resonant 
NWs, let us show first in Fig. [Tljthe classical plane wave 
refraction through a rectangular prism with an angle of 
6 = 30°, which actually involves two interfaces termi¬ 
nated in different photonic crystal directions (see the ap¬ 
pendix for details of the FEM numerical calculations). 
We can see how the wavefronts leave the prisms above 
the normal, a clear signature of a NIM medium. Indeed, 
the full simulation of the wave propagation (see Video 
1 in Ref. 02j Supplementary Information) revealed that 
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FIG. 10. (Color online) Top: Transmission through a slab 
made of a hexagonal array of Ag@Si core-shell nanowires (a as 
350 nm, Ri n = 90 nm, R out =170 nm), of the electromagnetic 
fields produced by in-plane point (horizontal/vertical) dipole 
sources at 235 THz (1.28 /im), as in Fig. [5j but for a slab 
surface along the TK direction. Bottom: Corresponding EFS 
(blue curves) of the hexagonal lattice along with that of light 
impinging from vacuum/air (black circle). A dashed green 
(vertical) line denotes wavevector conservation along the TK 
direction for a particular angle of incidence, producing ac¬ 
cordingly two refracted beams as indicated by phase/energy 
(blue/red, respectively) vectors. 


the wavefronts travel backward inside the prism (negative 
phase velocity), characteristic of a system having simul¬ 
taneously negative permittivity and permeability. The 
homogeneous prism with the retrieved effective parame¬ 
ters nicely reproduces the expected negative refraction, 
thus further supporting the accuracy of the effective pa¬ 



FIG. 11. (Color online) (a) Refraction through a 6 = 30° 
prism made of a NIM built with Ag@Si core-shell NWs with 
R in = 80 nm and R 0ut = 170 nm arranged in a hexago¬ 
nal lattice (a = 350 nm), illuminated by a transverse elec¬ 
tric (truncated) plane wave (1.215 /xm) parallel to the left 
interface: contour maps of the out-of-plane, z-component 
of the magnetic held amplitude, (b) Same as in (a), but 
for a homogeneous medium with e e ff = —0.45 + x0.04 and 
/jL e s = —1.7 — z0.21, so that n e s = —0.88 + z0.01 ( f.o.m . ~ 80) 
as obtained in Fig. [4] 


rameters and the NIM behavior of the NW array. 

We now introduce disorder in the prism geometry used 
in Fig. ED In essence, some of the constituents are either 
shifted or removed from their original position in the or¬ 
dered lattice. In Fig. [j~2^a) we explicitly show the NWs 
that have been shifted from their array positions in the 
prism and in Fig. |T2](b) the plane wave propagation into 
the disordered prism. Although the wavefronts are some¬ 
what perturbed by disorder upon crossing the prism and 
in the very near field region, we still observe the two key 
feature that characterize the negative index of refraction: 
an angle of transmission above the normal of the surface 
and a backward propagation of the wavefronts (negative 
phase velocity). 

Moreover, such NIM features are preserved even if 
stronger scattering disorder is introduced, as shown in 
the lower panel in Fig. |T2j where some NWs have been re¬ 
moved exactly along the wavefront propagation inside the 
prism. Actually, wave distortion is substantially larger 
than in the preceding case (shifted disorder), clearly ob¬ 
served in the near field close to the outgoing prism inter¬ 
face, although the backward propagation inside the prism 
and the (not so) far-field negatively refracted wavefronts 
are preserved. Bear in mind that, apart from positional 
disorder, if the NW filling fraction is substantially dimin¬ 
ished, a decrease in the (absolute) values of the negative 
effective parameters is expected^. 


Finally, disorder is also introduced in the slab geome¬ 
try used in Fig.[5j again by shifting some of the NWs from 
their original position in the ordered lattice, as shown in 
Fig. [13] Note that, to make more relevant the impact of 
disorder, the shifted NWs are located in the central part 
of the slab, through which most of the dipole transmitted 
power is concentrated. A moderate amount of disorder 
(top) does not disturb the flat-lens focusing. Only when 
the number of shifted NWs increases (see Fig. 13, bot- 
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FIG. 12. (Color online) Negative refraction through disor¬ 
dered rectangular prisms with an angle of 0 — 30° built with 
Ag@Si core-shell nanowires with Ri n = 80 nm and R ou t = 
170 nm arranged in a hexagonal lattice (a = 350 nm) illumi¬ 
nated by a (transverse electric) plane wave parallel to the left 
interface at 1.215 pm. Disorder is introduced in two different 
manners: (top left) some of the NWs (darker) are shifted by 
an average of 1.5 nm; (bottom left) some of the NWs are 
removed. Right column: Corresponding contour maps of the 
amplitude of the magnetic field (out-of-plane) component. 


tom), the transmitted pattern is slightly altered, though 
the transmitted focal spot is still reproduced. Vacancy 
disorder (not shown here) is extremely sensitive in this 
case to the vacant position for obvious reasons, leading to 
strong scattering if placed nearby the intermediate focus 
region inside the slab, as expected. Nevertheless, we have 
shown that the properties of our closely packed, doubly- 
resonant NWs are not only due to the periodicity of the 
lattice, but do survive even if substantial (either posi¬ 
tional or vacancy-like) disorder is introduced, bearing in 
mind that filling fraction should be preserved. 


V. CONCLUDING REMARKS 

In conclusion, we have carried out a full description of 
the photonic band structure of doubly-resonant core-shell 
metallo-dielectric nanowire arrays, in order to demon¬ 
strate their behavior as low-loss, isotropic bulk optical 
NIMs, as reported in Ref.42] from the retrieved effec¬ 
tive medium parameters. This must be done, as pointed 
out in Ref. |50] to rule out that their negative refrac¬ 
tion behavior is mainly a photonic-crystal diffraction ef¬ 
fect that depends on the specific geometry sample. The 
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FIG. 13. (Color online) Transmission through a slab made of 
a hexagonal array of Ag@Si core-shell nanowires (a = 350 nm, 
Rin = 90 nm, R out = 170 nm, TM surface direction), of the 
electromagnetic fields produced by in-plane (horizontal, cen¬ 
ter, and vertical, right) point dipole sources at 235 THz (1.28 
li m), as in Fig.[5](Poynting vector only), introducing disorder 
in the region below the dipoles by shifting the position of the 
NWs (darker ones) shown in the left column schemes: 16 (top 
row) and 37 (bottom row) NWs are randomly shifted. 


2D photonic band structure has been presented, iden¬ 
tifying the second photonic band with negative-index 
properties. The resulting EFSs in the first Brillouin 
zone are circular for most of this band, thus indicat¬ 
ing the full (2D) isotropic character. In addition, such 
isotropic NIM behavior is very clearly evidenced through 
numerical simulations in canonical examples for finite 
geometries: Namely, point dipole transmission through 
a slab and negative refraction through a prism. Both 
cases are compared to the same configurations, but for a 
slab/prism characterized by homogeneous dielectric per¬ 
mittivity and magnetic permeability exactly given by 
the effective parameter retrieval (from S-parameters), in 
very good agreement. This further confirms the effective 
(highly isotropic) bulk NIM behavior, in turn explicitly 
manifesting the extremely low losses and suitable (close 
to 1) impedance. Finally, disorder is introduced in such 
finite prisms to support the robustness of the NW struc¬ 
tures as NIMs, inferring that Mie resonances (rather than 
Bragg scattering) govern to a large extent their optical 
properties. 

In this manner, we have rigorously demonstrated 
(and improved) the highly isotropic and ultra low-loss 
(. f.o.m. ~ 200) character of Ag@Si nanowire arrays, ex¬ 
hibiting indeed robust flat lensing with subwavelength 
resolution and dipole position independence, in good 
agreement with the expected behavior for a negative- 
index homogeneous slab. Our results and procedure 
serve to fully assess such fascinating properties as optical 
NIMs of the doubly-resonant core-shell metallo-dielectric 
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nanowire arrays, paving the way towards realistic fabri¬ 
cation thanks to the tolerance to specific geometric con¬ 
siderations and inhomogeneities. 

Appendix: FEM numerical calculations 

The numerical calculations of photonic band struc¬ 
tures and wave propagation through slabs were carried 
out by means of full electromagnetic simulations using 
a Finite Element Method commercial software (COM- 
SOL Multiphysics 4.4). For photonic band structures, 
an eigenfrequency study is applied to the unit cell upon 
Bloch boundary conditions to reproduce an infinite sys¬ 
tem. The calculated photonic band direction is deter¬ 
mined by the periodicity of the fc-vector. 

In the case of the metameterial slabs, the simulated 
space consisted of an hexagonal array of of core-shell 
nanowires of height H ^6.7 fi m and length L « 35 
/im (22x100 nanowires), representing the slab, centered 
in a rectangle three times thicker than the slab and with 
the same length. In addition, all the physical system is 
surrounded by an additional layer of 2 /im that is set as 
a perfectly matched layer (PML) to absorb all the out¬ 
going radiation. The source of the electromagnetic field 
is simulated by an in-plane electric point dipole. For the 
homogeneous slab, the simulated space is the same as the 
metamaterial slab, but replacing the periodic array by a 
rectangle of the same H and L, and applying a transition 
boundary condition to its boundaries, setting its thick¬ 


ness to 1 nm and its optical properties to those of the 
surrounding medium (air). 

In all cases, silicon was assumed lossless in the fre¬ 
quency range studied, with a refractive index of nsi = 3.5 
and silver material constant were taken from experimen¬ 
tal values fromP3. The electromagnetic properties of 
the homogeneous medium were those determined by S- 
parameter retrieval procedure as in Refsl^MO, and the 
rest of the space was set to be air ( n a i r = 1). 

The meshing was done with the program built-in al¬ 
gorithm, which creates a triangular mesh. For photonic 
band structure calculations, the mesh maximum element 
size (MES) was set to 7.4 nm and the maximum element 
growth rate (MEGR) to 1.1. In addition, opposite edges 
of the unit cell were meshed identically by using the Copy 
Edge feature to ensure the symmetry of the system. For 
the transmission of the electromagnetic field radiated by 
a dipole, the MES was set to 100 nm and the MEGR to 
1.3. 
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